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2.  Objectives 

2. 1 .  Solve  Model  of  Contaminant  Transport  at  Air  Force  Bases 

This  problem  was  brought  to  our  attention  by  Marc  Jacobs.  At  his  suggestion  we 
contacted  Ben  Fitzpatrick.  TJie  mathematical  model  is  a  path  integral  with  an  integrand 
which  involved  the  solution  of  a  stochastic  differential  equation.  Rather  than  using*^our 
high  dimensional  integration  software  on  a  finite  dimensional  approximation  to  the  path 
integral,  we  decided  to  attack  the  numerical  solution  of  the  continuous  model  directly. 

Our  initial  objective  is  to  develop  algorithms  for  solving  a  model  problem  in  one 
space  and  one  time  dimension  and  to  embed  them  in  software.  We  will  compare  the  new 
algorithms  with  the  existing  Monte  Carlo  software.  If  the  results  are  promising,  we  will 
attack  the  higher  dimensional  transpon  equation. 

2.2.  Breaking  Intractability 

Continue  to  develop  and  apply  a  theory  on  breaking  the  computational 
intractability  of  continuous  problems  in  many  variables. 

2.3.  High  Dimensional  Integration 

In  panicular,  study  the  important  problem  of  high  dimensional  integration. 
Theoretical  results  we  obtained  earlier  suggest  that  low  discrepancy  methods  may  be 
superior  to  Monte  Carlo  for  at  least  some  applications  of  high  dimensional  integration.. 

2.4.  Software 


Build  and  test  software  running  on  a  heterogeneous  cluster  of  workstations  to  test 
the  theoretical  predictions  of  Section  2.3. 

2.5.  Implementation  Testing  of  Continuous  Problems 

Develop  a  theory  of  implementation  testing.  That  is,  given  a  mathematical  problem 
and  a  program  which  purports  to  compute  its  solution,  test  whether  the  program  does 
indeed  solve  the  problem. 


3.  Status  of  Effort 

ril  use  the  same  numbering  as  in  Section  2. 

3. 1 .  Solve  Model  of  Contaminant  Transport  at  Air  Force  Bases 

We  have  programmed  a  new  algorithm  for  solving  this  problem  and  compared  it  with 
the  standard  Monte  Carlo  approach.  The  new  algorithm  is  superior  with  respect  to  both  accuracy  and 
speed. 

3.2.  Breaking  intractability 

We've  made  good  progress  on  new  algorithms  for  the  numerical  computation  of  path  integrals. 
The  first  paper  has  been  published  in  the  physics  literature.  A  second  paper,  on  Feynman-Kac  type 
integrals,  is  in  preparation  A  new  algorithm  has  been  discovered  for  Feynman-Kac  integrals.  It  has  two 


major  advantages  over  Monte  Carlo:  It  provides  us  with  a  worst  case  guarantee  rather  than  only  a 
stochastic  guarantee  and  it  has  a  better  rate  of  convergence. 

3.2.1  Bellman  Equation 

We  have  staned  to  investigate  the  numerical  solution  of  the  Bellman  equation.  Again, 
we  believe  that  the  new  algorithm  will  be  superior  to  Monte  Carlo  in  providing  a  worst  case 
guarantee  with  a  better  rate  of  convergence. 

3.3.  High  Dimensional  Integration 

We  stated  in  our  last  progress  repon  that  “our  results  on  low  discrepancy  deterministic  methods 
have  stin-ed  considerable  interest  and  led  a  number  of  researchers  to  write  papers  trying  to  explain  our 
surprising  results.  We’ve  decided  to  tackle  this  question  ourselves.  A  paper  is  in  process  which  we 
believe  explains  not  only  our  success  in  very'  high  dimensional  integration  but  can  also^be  used  to  solve 
Bellman's  equation  and  other  important  applications. 

3.4.  Software  for  High  Dimensional  Integration 

We  have  continued  on  testing;  the  results  are  reponed  in  Section  4. 

3.5.  Implementation  Testing  of  Continuous  Problems 

This  work  continues  Jointly  with  David  Lee.  Lucent.  Three  papers  have  appeared. 

4.  .Accomplishments/New  Findings 

We  describe  our  accomplishments  in  two  areas. 

4.1.  Breaking  Intractability 

A  new  algorithm  has  been  created  for  Feynman-Kac  integrals.  It  has  two  major  advantages  over 
Monte  Carlo:  It  provides  us  with  a  worst  case  guarantee  rather  than  only  a  stochastic  one  and  it  has  a 
better  rate  of  convergence.  A  paper  describing  this  work  is  in  progress. 

4.2.  Software  for  High  Dimensional  Integration 


We  tested  the  Generalized  Faure  Sequence  on  a  360  degree  dimensional  model  problem  from 
mathematical  finance  and  found  that  we  beat  Monte  Carlo  by  a  wide  margin.  We  and  others  believed  that 
this  was  due  to  the  non-i.sotropic  nature  of  the  model  problem. 

We  decided  to  test  a  hard  isotropic  problem  that  had  been  published  in  “Computers  in  Physics”  in 
1 996.  We  again  found  the  Generalized  Faure  Sequence  beat  Monte  Carlo  by  orders  of  magnitude.  Our 
results  appeared  in  “Computers  in  Physics”  in  November,  1997. 

These  results  are  empirical.  We  need  a  theory  which  tells  us  for  which  class  of  integrals  low 
discrepancy  methods  beat  Monte  Carlo.  We  believe  we  have  found  the  explanation  for  non-isotropic 
integrals;  a  paper  is  in  progress.  The  explanation  for  the  new  empirical  results  on  isotropic  integrals  is 
open. 
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Traub  -  FANNIE  MAE,  Washington,  DC,  October,  1995 

Traub  -  Courant  Institute  of  Mathematical  Sciences,  New  York  University,  New  York,  NY,  November. 


Traub  -  Computer  Science  Department,  University  of  Toronto,  Toronto,  Canada, 

November,  1995 

Traub  -  Computer  Science  Department,  University  of  Waterloo,  Ontario,  Canada, 

November,  1995 

Traub  -  Fields  Institute,  Toronto,  Canada,  November,  1995 

Traub  -  Address  to  Chief  Scientists  and  technical  Officers.  U.S.A.F.,  Santa  Fe.  NM.  July,  1996 

Traub  -  Workshop  on  Fundamental  Limits  to  Knowledge  in  Economics,  Santa  Fe 
Institute,  Santa  Fe,  NM,  August,  1996 

Traub  -  .Mter  dinner  address.  Century  Association.  New  York,  NY,  October,  1996 

Wozniakowski  -  Foundation  of  Computational  Mathematics,  co-organizer  and  speaker  of  the  workshop, 
“Information-based  Complexity”,  Rio  de  Janeiro,  January,  1997 

Wozniakowski  -  City  University  of  Hong  Kong,  Hong  Kong,  January,  1997 

Wozniakowski  -  Workshop,  “Computability  and  Complexity  in  Analysis",  Dagstuhl,  Germany,  .\pril. 


Wozniakowski  -  International  Computer  Science  Institute,  Berkeley,  CA,  July,  1997 
Wozniakowski,  Los  Alamos  National  Laboratory.  Los  .Alamos,  NM,  August,  1997 
7B  Advisory  Functions  to  Air  Force,  DOD,  and  Other  Laboratories 

7B.  1.  On  July  24,  1996  I  was  invited  to  attend  a  meeting  of  the  Air  Force  Chief  Scientists  Group  in 
Santa  Fe.  The  group  consisted  of  chief  scientists  and  chief  technical  officers  from  many  Air  Force  Bases 
and  Laboratories,  The  meeting  was  organized  by  Edwafd  Feigenbaum,  Chief  Scientist  of  the  Air  Force. 
Participants  included  Jo.seph  Janni  and  Helmut  Hellwig. 

7B.2.  On  August  12,  1996  I  participated  in  a  meeting  of  the  Highlands  Group  in  Santa  Fe  at  the 
invitation  of  Frank  B.  Horton.  Ill,  Principal  Deputy  Assistant  Secretary  of  Defense  (Command,  Control, 
Communications,  and  Intelligence)  and  Andrew  W.  Marshall,  Director,  Net  Assessment,  Office  of  the 
Secretary  of  Defense.  We  worked  on  scenarios  regarding  the  impact  of  information  technologies  on  the 
military  for  a  time  horizon  of  up  to  twenty  years. 

7B.3.  In  November,  1996  I  participated  at  a  meeting  of  the  Highlands  Forum  at  the  Naval  Academy, 
Annapolis,  to  di.scuss  information  warfare.  We  worked  with  members  of  the  Presidential  Commission  on 
Protecting  the  National  Information  Infrastructure. 

7C.  Transitions 

We  have  assigned  the  rights  to  license  our  high  dimensional  integration  software 


system  FINDER  to  Columbia  University.  Columbia  has  licensed  the  software  to  FANNIE 
MAE.  The  University  is  actively  seeking  other  licensees. 

8.  New  Discoveries 

We  have  applied  for  a  patent  on  our  software. 
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